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ABSTRACT 

A simple, yet powerful, algorithm for computed tomography of the solar corona is presented and demon- 
strated using synthetic EUV data. A minimum of three perspectives are required. These may be obtained from 
STEREOfEUVl plus an instrument near Earth, e.g. TRACE or SOHOfETT. 
Subject headings: methods: data analysis — Sun: corona — techniques: image processing 



1. INTRODUCTION 

Observations by the Extreme Ultraviolet Imager (EUVI, 
Wuelser et al. 2004) aboard STEREO provide the first si- 
multaneous, stereoscopic image pairs of the solar corona and 
transition region. Ideally, these data are simple projections 
through an optically thin corona. However, the 3D distribu- 
tion of emission is difficult to estimate with only two pro- 
jections (see Gary et al. 1998, and references therein). This 
difficulty may be explained as follows. Consider an object 
I(x,y,z) on domain D. The coordinates x and y need not be 
orthogonal, but z is orthogonal to x and y. Given two projec- 
tions f(x,z) — fpldy and g(y,z) — \ D ldx, a plausible recon- 
struction of / is 



f(x,z)g(y,z) 
T(z) ' 



(1) 



where T(z) — f D I(x,y)dxdy = J D fdx = J D gdy is the total 
emission of a plane of constant z- Unfortunately, this solution 
fails utterly in practice. For each pair of sources in /, /' intro- 
duces a pair of "ghost" artifacts. These are systematic errors, 
independent of the noise and apparently unavoidable. Tradi- 
tional regularization strategies are not fruitful: /' is positive, 
is as smooth as the observed images / and g, and is precisely 
the maximum entropy solution. More information is therefore 
required to guide the tomographic reconstruction. 

Previously described approaches to the STEREO coronal 
tomography problem rely on assumptions about the geom- 
etry of the coronal plasma distribution. The triangulation 
method (Gary et al. 1998; Aschwanden 2005; Aschwanden 
et al. 2008) assumes loops with circular cross-section, and re- 
lies on the identification of the same loops in both images. 
The magnetic tomography approach (Wiegelmann & Inhester 
2006) also assumes loops with circular cross-section, and in- 
corporates magnetic field extrapolations to constrain loop ge- 
ometries. These methods are powerful, but they require as- 
sumptions about things that one might reasonably hope to 
learn from the 3D reconstruction. 

I propose that EUV images taken from a third perspective — 
e.g. TRACE or SOHO/EYT — may provide adequate additional 
constraints for coronal tomography, without any assumptions 
about loop geometry or magnetic fields. I describe a simple 
computed tomography algorithm, fast enough to run in real 
time, and demonstrate its performance using synthetic data 
with three viewpoints. 

2. ALGORITHM 
Electronic address: kankel@solar.physics.montana.edu 



The Smooth Multiplicative Algebraic Reconstruction Tech- 
nique (SMART) presented here has been developed to solve 
a mathematically analogous problem of reconstructing spec- 
tra for the MOSES sounding rocket payload (Kankelborg & 
Thomas 2001; Fox et al. 2003). Iterative multiplicative alge- 
braic reconstruction techniques (MART), perhaps inspired by 
the separable solution to the two-view problem (equation 1), 
have been available for many years (Okamoto & Yamaguchi 
1991; Verhoeven 1993). Gary et al. (1998) used MART along 
with volume constraints derived from magnetic field extrapo- 
lations to reconstruct a pair of loops from two simulated XUV 
images. The unique algorithmic features of SMART are iter- 
ative smoothing and an adaptive correction strategy. These 
refinements improve numerical stability and promote conver- 
gence to a compromise between smoothness and goodness- 
of-fit, leading to a reduced chi-squared of unity. 

In the A^-view tomography problem, an object I(x,y,z) on 
domain D is known only by N projections f m , taken at angles 

a .1 
Cm- 



f m (x,z) = / M(9 m )I(x,y,z)dy. 



(2) 



The operator &(0 m ) rotates the object / by angle 9 m about the 
z axis. 2 SMART uses the projections f m to estimate I(x,y,z) 
by the following steps: 

1. Create an initial guess, I'(x,y,z) = 1 on D. 

2. Initialize correction weights, y m = i . 

3. /' <— /' *K (smoothing kernel K defined by eq. 7). 

4. Calculate projections fj„(x,z) = f D &(9 m )I' dy. 

5. Calculate correction factors, 



C m (x,y,z)=£(-6 m ) 



fm(x,z)' 



(3) 



Note that a nontrivial y-dependence is introduced 
through the rotation. 



6. Apply corrections weighted by y m 



(4) 



1 In our coordinates, 0,„ is a light-handed rotation of the object; it therefore 
corresponds to the eastward heliographic longitude of the observer. 

2 Rotation S could incorporate compound angles with altitude, azimuth 
and roll. The extension of SMART to the general case is straightforward. 
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7. Calculate reduced chi-squared for each projection, 



%R,n 



1 



(/-/') 



>\2 



N x N y N z ~ z 



8. Adjust correction weights, y m , using a control law de- 
signed to drive %\ m — » 1. 

9. Repeat steps 3-8 until converged. 

The /' array is oversized so that rotations will not move 
any of the emission outside the volume. The subphotospheric 
portion of the volume is zeroed in step 1, and will remain zero 
since the corrections are multiplicative (eq. [4]). 

Since the correction factors (eq. [3]) are ratios of positive 
numbers, the reconstruction is always positive. The strength 
of the mth applied correction is governed by y m . Step 8 sets up 
feedback to establish a dynamic equilibrium between smooth- 
ing and correction, so that %\ m tends toward unity. Our control 
law, which has two adjustable parameters (a,b), modifies y m 
for iteration n + 1 using a linear combination of the previous 
and current values of j| m : 



yi+i 

hn 



X,Rm,n 



(5) 
(6) 



The algorithm is implemented in IDL, with rotations &(d m ) 
perfomed via cubic convolutional interpolation. It is possible 
that spurious negative voxels could be introduced during the 
rotation, but negative values are eliminated from our projec- 
tions by thresholding: f m > 0. 

The normalized smoothing kernel, K(x,y,z), is defined on 
the discrete space of voxels as follows: 



K ijk - 



„(i+j+k) 



l+s(l+4c+12c 2 + 8c 3 



(7) 



This form of K is not crucial, but it is designed to be very 
nearly isotropic. The adjustable smoothing parameter, s G 
(0, 1], affects the rate of convergence but has no discernible 
effect on the result. 

3. SYNTHETIC DATA 

A volume of synthetic coronal emission was prepared as 
a test target for the SMART algorithm. The test target is 
more complex than any that have proven tractable for pre- 
vious approaches to tomographic analysis for STEREO. It is 
not an attempt at detailed atmospheric modeling, but resem- 
bles a small active region. I began with a potential field de- 
fined by four sub-photospheric magnetic charges. The result- 
ing line-of-sight photospheric magnetic field is shown in fig- 
ure 1. The field lines in the figure emerge from five small, 
square "heated" patches in the postive polarities on the pho- 
tospheric plane. A large number of field linese were traced 
from random points in the volume to both footpoints (or to 
the boundary). Emission was added to each voxel crossed by 
the field line, in proportion to the arbitrary heating value at its 
positive photospheric footpoint. Figure 2 shows the volume 
model projected along each of the three axes of figure 1 . The 
coronal flux tubes have complicated geometries, including a 
broad, vertical fan of emission above the postive pole that is 
associated with a coronal magnetic null at a height of ~ 53 
grid points. All images of EUV emission in this letter are 
square-root scaled to bring out faint features. 




FIG. 1. — Vertical magnetic field (shaded plane) and field lines within the 
emission from the model corona, five "Heated" patches within the positive 
photospheric magnetic poles are colored in black. Note the sharp deflection 
of the upppermost field line from the coronal null. 
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FIG. 2. — The model coronal volume viewed along each of the three axes 
of figure 1 . Intensities are square-root scaled. 

SMART was tested on numerous synthetic observations of 
the model coronal volume, each time placing the model ac- 
tive region at a different random northern heliographic lat- 
itude over [0°,30°], a random heliographic longitude over 
[-40°, 40°], and a random tilt over the interval [-180°, 180°]. 
For simplicity, in our coordinate system the solar equator co- 
incides with the equatorial plane. Observations were pro- 
jected for three distant virtual instruments in the ecliptic, ob- 
serving from heliographic longitudes — 40°, 0°, and 40°. The 
twin STEREO spacecraft will reach similar separation an- 
gles in October, 2008. The images were normalized so that 
the brightest pixel among the three images had 3000 counts. 
The 9 = (-40°,0°,40 o ) projections shown in figure 3 cor- 
respond to an example observation with the region placed at 
latitude 16.51° N, longitude 14.65° W, and tilt 29.63° counter- 
clockwise. Poisson noise was applied to the images prior to 
passing them to the SMART algorithm for inversion. Inten- 
sities are square-root scaled to show the noise more clearly. 
The mean value of the nonzero pixels is 216 counts. 



Coronal Tomography 



3 



13 = -40° 13 = +00° 13 = +40° 



FIG. 3. — Synthetic data taken from three virtual spacecraft. Intensities are square-root scaled. 



4. INVERSION RESULTS 

The synthetic data were inverted using 15 SMART itera- 
tions. This was typically sufficient to converge to Xr m = 
1 ±0.01 . The rate of convergence is affected by the smoothing 
parameter and by the two adjustable parameters in the con- 
trol law. The examples shown in this Letter used s = 0.5, 
a = 0.05, b = 0. 16. Within broad limits, the results are insen- 
sitive to these parameters. For example, if the smoothing is 
removed altogether, numerical instabilities arise; if it is made 
too strong, then it will not be possible to reach %\ m = !• F° r 
our computational volume of 169 3 voxels, an inversion runs 
in a few minutes on a laptop computer. 

The simulated data in figure 3 gives rise to an accurate re- 
construction. 

Figure 4 shows reconstruction results for the data in figure 
3, compared to noise-free visualizations of the coronal vol- 
ume model. The Earth viewpoint at the top of figure 4 cor- 
responds to the middle panel of figure 3. Comparing these 
two figures shows that the photon noise has been largely sup- 
pressed in the SMART reconstruction. The success at re- 
covering 3D geometry is best illustrated by the east-west and 
south-north projections, which are 50° and 90°, respectively, 
from the nearest available observing angles in the synthetic 
data. All of the true features have been recovered. Square- 
root scaling helps to bring out the artifacts, which are few and 
faint. There is slight blurring, and minimal ghosting. These 
results are typical of hundreds of realizations tried so far. 

A second example, with loops nearly parallel to the ecliptic 
plane, is given in figure 5. The horizontal loop orientation is 
very challenging because only the ends of horizontal features 
provide any depth cues. Altough views from within the eclip- 
tic plane are reproduced well, the example shows relatively 
poor reconstruction of an out-of-ecliptic view (lower panel). 
Animated versions of figures 4 and 5 are provided in the elec- 
tronic version of the Journal. 

5. DISCUSSION AND CONCLUSIONS 

Coronal tomography is possible with as few as three view- 
points, making no prior assumptions about coronal morphol- 
ogy or magnetic fields. The test cases demonstrate recovery of 
complex geometry without reference to magnetic field extrap- 
olations or assumptions about loop geometry. Loops that run 
in an east-west direction, however, provide insufficient depth 
cues for three instruments confined to the ecliptic. 

The SMART algorithm described here provides a noise- 



insensitive tomographic reconstruction by finding an opti- 
mal balance between goodness of fit and local smoothness. 
STEREO will obtain data at large separation angles in Fall 
2008. Observations from SOHO and/or TRACE at that time 
should allow the best opportunity for application of SMART 
to coronal tomography. 
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True SMART 
Earth view 
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From east, looking west 



From south, looking north 



FIG . 4 . — Comparison of true and reconstructed coronal volumes from three 
orthogonal points of view. Intensities are square-root scaled. 
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From south, looking north 



FIG. 5. — Same as figure 4, but with loops nearly parallel to the ecliptic. 
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